Neutrino mixing and mass hierarchy in Gaussian landscapes 
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The flavor structure of the Standard Model may arise from random selection on a landscape. 
In a class of simple models, called "Gaussian landscapes," Yukawa couplings derive from overlap 
integrals of Gaussian zero-mode wavefunctions on an extra-dimensional space. Statistics of vacua 
are generated by scanning the peak positions of these wavefunctions, giving probability distributions 
for all flavor observables. Gaussian landscapes can account for all of the major features of flavor, 
including both the small electroweak mixing in the quark sector and the large mixing observed in 
the lepton sector. We find that large lepton mixing stems directly from lepton doublets having 
broad wavefunctions on the internal manifold. Assuming the seesaw mechanism, we find the mass 
hierarchy among neutrinos is sensitive to the number of right-handed neutrinos, and can provide a 
good fit to neutrino oscillation measurements. 



I. INTRODUCTION 

The Standard Model of particle physics, taken to in- 
clude neutrino masses, is described by a quantum field 
theory with about two dozen input parameters. The fla- 
vor observables — quark and lepton masses, mixings and 
CP phases — constitute a large fraction of these inputs, 
reflecting the fact that we have not found a simple princi- 
ple (like gauge coupling uniflcation) to relate all of them. 
It is possible that there is no deep meaning behind the 
precise values of many of the flavor parameters that we 
observe. Furthermore, the set of Yukawa couplings in our 
vacuum may be one among many possibilities that the 
fundamental (microscopic) theory admits. If this is the 
case, then one can understand the measured values of fla- 
vor observables only in terms of statistical distributions; 
still the large number of flavor observables gives hope for 
a discerning statistical analysis. These ideas have been 
pursued in Refs. Note that in string theory, 

flux compactiflcation succeeds in stabilizing some mod- 
uli, the statistics of fluxes in the internal space generates 
statistics of moduli values, and these may in turn gener- 
ate statistics of Yukawa couplings Thus, the current 
understanding of string theory supports a statistical pic- 
ture of flavor. 

The Yukawa couplings of the Standard Model can- 
not be completely random numbers. The eigenvalues of 
Yukawa matrices are hierarchically separated, at least in 
the quark and charged lepton sectors, and the charged 
electroweak current pairs each left-handed up- type quark 
almost uniquely to a left-handed down-type quark; fur- 
thermore this pairing combines the heaviest up-type 
quark to the heaviest down-type quark, and similarly for 
the middle and lightest quarks. These patterns are re- 
spectively referred to as hierarchy, pairing structure, and 
generation structure, and a statistical theory of flavor 
observables should explain each of them. 



An intuitive explanation of these patterns is available 
in string theory compactiflcation Fields in the low- 
energy effective theory (e.g. the Standard Model) are 
Kaluza-Klein zero modes on some internal manifold, and 
Yukawa couplings are calculated by overlap integration of 
the zero-mode wavefunctions of the three flelds relevant 
to any coupling. Overlap integration of localized wave- 
functions generates hierarchy [6|. Localized wavefunc- 
tions of the quark doublets and the Higgs boson introduce 
correlation between the up-type and down-type Yukawa 
matrices, giving rise to pairing structure and generation 
structure. Meanwhile, localized zero-mode wavefunctions 
are fairly easy to obtain in torus-flbered compactifica- 
tions with small torus fiber, see e.g. Refs. These basic 
features of quark flavor are nicely explained in Gaussian 
landscapes [1, [1] , which are proposed as toy models of 
the landscape, capturing the essential features of torus- 
fibered compactification of Heterotic string theory and 
its dual descriptions. 

The observation of large mixing angle neutrino oscilla- 
tions has been welcomed with surprise, as it reveals that 
the matter fields of the Standard Model are not simply 
three copies of a spinor representation of SO (10). The 
lack of pairing structure in the lepton sector presumably 
indicates that flavor structure is quite different between 
quark and lepton sectors. How do the quark and lep- 
ton sectors come to have different flavor structure, and 
why? These questions point to a deeper, more funda- 
mental understanding of the microscopic origin of flavor. 
Such an understanding might influence the way we study 
leptogenesis, or how we make predictions for the flavor 
observables (such as 0™^^ and CP-violating phases) to 
be measured in future neutrino experiments. 

Conventional theories of flavor assume that mass hier- 
archies and mixing angles are determined by approximate 
flavor symmetries. With the discovery of large neutrino 
mixing angles, and at most a modest neutrino mass hier- 
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archy, Ref. [l| proposed a statistical "neutrino anarchy" 
in the lepton sector, with lepton doublets not determined 
by any symmetry. In the Gaussian landscape [3, a 
microscopic origin for quark mass hierarchies and mixing 
angles emerges from the localization of zero modes on 
an internal manifold. It was a natural guess in Ref. Q 
that large lepton mixing angles would follow from lepton 
doublets having less localized zero-mode wavcfunctions 
than quark doublets, allowing a description of "/lo?/;" the 
quark and lepton sectors differ. 

We here conclude that this guess is indeed correct. The 
present authors claimed in Refs. 0,01 that large mixing 
in the lepton sector required both non-localized wave- 
functions of lepton doublets and large complex phases. 
The latter were seen as necessary to prevent a tidy can- 
cellation of terms in the PMNS matrix. However, we 
have since found that the above cancellation proceeds 
as an artifact of ignoring kinetic mixing of fermions in 
Refs. With kinetic mixings taken into account in 

Gaussian landscapes, we find that complex phases are 
not necessary for large lepton mixing. 

We also find that the (statistical distributions of) low- 
energy neutrino masses, generated by the see-saw mecha- 
nism in Gaussian landscapes, are sensitive to the number 
of right-handed neutrinos. When there is a large num- 
ber of right-handed neutrinos, their mass eigenvalues be- 
come more densely packed, and the hierarchy between 
the Am^ of solar and that of atmospheric neutrino oscil- 
lations need not be large. The existence of many right- 
handed neutrinos is a natural consequence of string the- 
ory compactification, because right-handed neutrinos are 
just SU(5)GUT-singlet moduli fields, and there are often 
many moduli. Note that a dense concentration of right- 
handed neutrino masses also implies that more than one 
neutrino might contribute to thermal leptogenesis. 

The remainder of this paper is organized as follows. 
The Gaussian landscape is set up in Section [TTl where we 
extend the construction of Refs. to account for over- 
lap integrals that may lead to non-canonical kinetic terms 
in the low-energy theory. Our main results appear in Sec- 
tion mil where we describe a numerical simulation of a 
Gaussian landscape and comment on its (in)sensitivity to 
model assumptions. In Section [TVl we briefly outline how 
certain analytical approximation methods, developed in 
Ref. [3l to provide insight into a Gaussian landscapes, 
can be extended to our revised construction. Concluding 
remarks are given in Section |Vl 



II. GAUSSIAN LANDSCAPES 

A. Kaluza— Klein Reduction 

Consider a low-energy effective theory (e.g. the Stan- 
dard Model) obtained by compactiiying a supersymmet- 
ric Yang-Mills theory of group G D SU(3)c x SU(2)l x 
U(1)f = GsM on a higher-dimensional spacctime. Let 
the internal manifold have D dimensions, and denote 



it by X. Wc assume a gauge- field background exists 
on X such that the unbroken gauge symmetry is re- 
duced from G to Gsm- Fermions in the low-energy ef- 
fective theory may be identified with the Kaluza-Klein 
zero modes of gauginos. The irreducible decomposition 
of Q under Gsm may contain (3, 2)+^/*^ -I- (3,l)-2/3 + 
(3,l)+i/3+ (l,2)-i/2+ (l,l)+i of Gsm, i.e. a set of 
representations corresponding to a "generation" of Stan- 
dard Model fermions. In each irreducible component, the 
higher-dimensional gaugino '^{x,y) has a Kaluza-Klein 
decomposition 

vi/"(x, y) = ^ C(^) v'^iy) + E ^?(^) ' (1) 

i I 

where a labels the irreducible representations listed 
above (a ~ q, u, d, £, or e), i runs over the Kaluza- 
Klein zero modes, and / runs over the massive modes. 
The wavefunctions 'p{y) denote mode functions of the 
Kaluza-Klein decomposition, while the ''p{x) correspond 
to fields in the 3+1 dimensional effective theory below the 
Kaluza-Klein scale [x denotes coordinates in 34-1 dimen- 
sional Minkowski space, y denotes those of the internal 
space X). Note that the i^ix) arc spinors of S0(3,l), 
while the ^(y) are spinors of SO(D), but spinor indices 
are suppressed. 

If the topology of X and the gauge-field background on 
X are chosen appropriately, then there are three Kaluza- 
Klein zero modes in each irreducible representation, just 
like in the Standard Model. In what follows, we only 
consider vacua with such properties. If the irreducible 
decomposition of g contains a singlet of Gsm, then zero 
modes of this component may be identified with right- 
handed neutrinos. Wc extend the range of the label a to 
include these right-handed neutrinos. Because we do not 
require an unbroken SO(IO) symmetry in the low-energy 
effective theory, the number of right-handed neutrinos is 
not necessarily three. In fact, the right-handed neutri- 
nos are supersymmetric partners of gauge-field moduli in 
compactifications with supersymmetry, in which case it 
is quite likely that there are many of them. 

The kinetic terms of the three fermions in each repre- 
sentation arise from dimensional reduction of the gaugino 
kinetic term in higher dimensions: 

where 

K?j = M^J^ d^y^g Tptiv) V-iv) ■ (3) 

Here A/* is the cutoff of the 4+D dimensional super 
Yang-Mills theory, and again D is the number of extra 
dimensions. 

Meanwhile, the Higgs doublet of the Standard Model 
may be identified with the Kaluza-Klein zero mode of 
a gauge field Am{x,y) of g on X. The Kaluza-Klein 
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decomposition of Am{x, y) may take the form 

A^{x, y) = h{x) ^tiy) + '*^(^) + l^-c- (4) 

I 

Here h[x) is a complex scalar field of the 3+1 dimen- 
sional effective theory and is in the (1, 2)+^/^ repre- 
sentation of GsM- The zero-mode wavefunction of this 
scalar is (p'^(y). The Hermitian conjugate part contains 
h*{x) (pjjj (y)*. The Yukawa couplings of quarks and lep- 
tons then arise from dimensional reduction of the super 
Yang-Mills interactions on X: 

Jx 

^A^''C(^)M2;)V','(^), (5) 

where 

K' = 9*M^I d''y^^'i{y)T^^t{y)^]{y), (6) 
Jx 

for the up-type and neutrino Yukawa couplings, while 
replacing h[x) with h*{x) in Eq. dSl) and 95j'„(y) with 
'Pmiu)* Eq. ([6]) for the down- type and charged-lepton 
Yukawa couplings. Here is a dimensionlcss cou- 
pling constant of the gauge interaction of G on the 
(£'+4) dimensional spacetime. This is essentially how 
the Standard-Model Lagrangian is obtained in the com- 
pactification of Heterotic / Type I string theories. Fur- 
thermore, some vacua of Type HA / M-theory / Type 
HB / F-theory compactifications are also approximated 
(albeit poorly) by this description, because of string du- 
ality. 

The Kaluza-Klein zero modes of any irreducible com- 
ponent of the gaugino form a vector space, because the 
massless Dirac equation on a D dimensional manifold 
X is linear in (p{y). One is free to choose any basis 
for the vector space of these zero modes. The kinetic 
mixing coefficients K°j and the Yukawa couplings Xfj in 
the low-energy effective Lagrangian look different for dif- 
ferent choices of basis, but this is simply due to field 
redefinitions — observables are independent the choice of 
basis. Thus, when an ensemble of {Kfj, Xfj} is deter- 
mined by considering vacua with various geometries and 
gauge-field configurations, the statistics of the flavor ob- 
servables are basis independent. 



the generic possibilities. From a "bottom-up" perspec- 
tive, one might first identify our corner of the landscape 
using phenomenological considerations, and only after- 
wards try to explain why this subset of vacua is preferred, 
whether it be by sheer statistics, dynamical effects, or an- 
thropic selection. While we do not attempt to identify 
all of the subsets of vacua in which our vacuum is typi- 
cal, the present authors pointed out in Refs. [1, 3 that 
the observed patterns of flavor are typical in at least one 
particular corner of the landscape. 

If X is a T^-fibration on a {D — m) dimensional base 
manifold B, and if the typical size of T*" is smaller than 
that of i3, then it is known that a Kaluza-Klein zero 
mode is effectively localized in a (Z) — 2m) dimensional 
subspace of B. This localization can be understood in- 
tuitively as follows: gauge fields tangential to the re- 
directions effectively become m independent mass pa- 
rameters for fields on B, and if the values of these "mass 
parameters" vary over _B, then fermion zero modes are lo- 
calized on B where the mass parameters vanish — a mech- 
anism known as a domain wall fermion. Since there are 
TO independent mass parameters, the zero-mode fermions 
are localized in m directions within the (D — m) dimen- 
sional manifold B. Wavefunction amplitudes decrease 
rapidly away from the locus of localization; the wave- 
functions are approximately Gaussian in profile.^ 

As pointed out in Ref. localized wavefunctions 
on extra dimensions easily generate hierarchically small 
Yukawa eigenvalues. Refs. [1,13] found this result to read- 
ily translate to the compactification scenario described 
above. Furthermore, it was found (and will be elaborated 
upon below) that this picture explains the other major 
patterns of flavor, including the pairing structure and 
generation structure in the quark sector, and large lep- 
ton mixing. Thus, "Gaussian landscapes" are proposed 
as toy models of our corner of the string landscape. Un- 
derstanding Gaussian landscapes can shed light on the 
microscopic dynamics relevant to our corner of the land- 
scape, and thereby help address the question "wftj/" of 
the introduction. 



C. Gaussian Landscapes 

In practice, it is extremely technically involved to carry 
out the program outlined in Section [II Al One must ffi'st 



B. Gaussian Zero- mode Wavefunctions 

If one considers generic internal manifolds X and 
generic configurations of gauge-field backgrounds on the 
various X, then one expects generic Yukawa couplings 
to arise in the low-energy effective theory. Given that 
we observe patterns in the flavor observables (such as 
those listed in the introduction), we expect our cor- 
ner of the landscape to be described by some subset of 



^ There is a phenomenological motivation to consider the local- 
ized gauge-field flux of U(l)y ■ When such a U(l)y flux is 
localized in the internal space X, the zero-mode wavefunctions 
tend to have (often linear-exponential) hypercharge-dependent 
special behavior around the localized flux. As a toy model, the 
Gaussian landscape below does not include this effect. This is in 
part because the localized U(l)y flux is just one of many pos- 
sible ways to break an SU(5)gut symmetry. Although it may 
be possible to modify the Gaussian landscape to implement this 
mechanism, we do not do this here. 
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identify a stable gauge-field configuration on a compact- 
ified space X. This involves solving a non-linear partial 
differential equation on a curved manifold. One must 
then find the zero mode wavefunctions Lp°'{y); note that 
(p°'{y) is a multi-component field, because it sits in a non- 
trivial representation of SO{D) as well as in a non-trivial 
gauge-field background. Finally, an overlap integration 
must be performed over X to obtain K^^ and Xf^. Note 
that one must calculate the metric on X in order to per- 
form this overlap integration, because the metric enters 
Eqs. (131) and ©. AH this generates just one element of 
the landscape ensemble of {K°-j, Xfj}; the process must 
be repeated many times to get a sense of the statistics of 
the landscape. Ref. Q takes a step toward performing 
such an analysis. 

However, all of this may not be necessary to under- 
stand the observed patterns of flavor. If we may speak 
with the benefit of hindsight, the distributions of fiavor 
obscrvables are quite broad across subsets of the land- 
scape ensemble. Our vacuum is selected randomly (mod- 
ulo any anthropic or cosmological effects) from the land- 
scape, so it is not important to understand observables 
with great precision. Instead, we seek to understand the 
broad-brush patterns of fiavor, assuming that the precise 
values of flavor observables are essentially accidents. 

For this purpose, we introduce a number of simplifying 
assumptions, hoping to capture the essence of flavor in a 
certain corner of the string landscape, while setting aside 
the complicating details. First, we use only the base man- 
ifold B when calculating overlap integrals to construct 
JsT?- and Afj*. Second, although the zero- mode wavefunc- 
tions if"" are multi-component fields, we treat them as if 
they were single component fields. Since we consider the 
approximately Gaussian profiles of zero-mode wavefunc- 
tions to be the origin of hierarchical Yukawa couplings, 
we expect single component wavefunctions to carry the 
most important information. For concreteness, we use 
wavefunctions with strictly Gaussian form (made peri- 
odic on B); e.g. for the ith zero mode of the representa- 
tion a. 



ipfiy) ~ exp 



1 + irg 
2dl 



(7) 



The peak position is a point in B, the width da deter- 
mines the degree of localization of the Gaussian profile, 
and Ta is a possible complex phase that we include for 
later reference.^ The normalization of these wavefunc- 
tions is arbitrary (it does not affect observables). 



^ The inclusion of tlic complex phase ra is not without motivation; 
it is an approximate form for zero-mode wavefunctions on a torus 
with generic complex structure (see Section 7 of Ref. for more 
details). However, it is unclear how complex phases enter the 
zero-mode wavefunctions when B is not a torus. We introduce 
complex phases into zero-mode wavefunctions in this way merely 
to explore a possibility. 



To establish a more convenient notation for later, we 
write the relevant terms of the low-energy Lagrangian, 

C D aiij^'Df, aj + A-'^ Ui h + Xfj di qj h'^ 

+ Ay- e^ Ij ft-t + A-j Ui Ij h + Af/ n^ fij , (8) 



where ni denotes a (heavy) right-handed neutrino. The 
indices i,j label generation, but note that we do not as- 
sume three "generations" of right-handed neutrinos. The 
kinetic matrices are given by overlap integration of 
zero- modes, c.f. Eq. here written 



K^^Mt" d^v^g^r{y)v';{y)- 

J B 

The up-type Yukawa matrix, c.f. Eq. is given by 



(9) 



Ay 



5.A/f / d^yV~9P{y)^ny)^]{y)v\y), (lo) 



where for later reference we have introduced a phase func- 
tion p{y) that might arise due to the "F-structure" in 
T"^^'^ of Eq. ©.3 The Yukawa matrices A^^ , A^^- , and A^; 
are all defined analogously. Although the wavefunction 
of the Higgs boson will in general have complicated struc- 
ture, in the Gaussian landscape we treat it as a single- 
component wavefunction with Gaussian profile.^ This is 
reasonable, since linear changes in the gauge-field back- 
ground result in quadratic changes in the mass-square 
parameter in the quadratic mode equation of the bosonic 
zero modes; note that the ground-state wavefunction of 
a harmonic oscillator is Gaussian. 

We assume low-energy left-handed neutrino masses 
arise from the seesaw mechanism. Hence these masses 



^ To be more precise, we can write F"" = T^e^, where the 
are fixed- value matrices, like the Pauli matrices, that satisfy the 
algebra of the internal SO{D). Just like with the Pauli matri- 
ces, the PP may have mutually different complex phases, e™ 
is a vierbein on the internal manifold, and thus depends on y. 
Therefore, we expect a y-dependent complex phase in P'". Our 
introduction of a phase function p{y) is an attempt at modeling 
the possible effects of this. 

* Note, however, that there may be important consequences of the 
wavefunction ip^{y) (and </'m (?/)*) having multiple components. 
Suppose that G D G' X Gsm^ a-nd that the G' symmetry is bro- 
ken by the gauge-field background. Then the ^1{y) are all zero 
modes under the gauge-field background in G', each in a certain 
representation pa- Suppose also that ip^{y) is in a representa- 
tion ph, and ¥'m(?/)* in . Then the up- type Yukawa couplings 
are given by overlap integration of trivial components in the irre- 
ducible decomposition of pa=u ® Pa—q ® Ph > while the down-type 
Yukawa couplings by that of trivial components of Pd^ Pq Ph ■ 
Thus, a single component wavefunction ip^{y) for A" and A" and 
precisely its complex conjugate wavefunction ip'^{y)* for A'' and 
X'^ may introduce too restrictive a relation between the Yukawa 
couplings in the different sectors. To account for this issue, one 
may extend the Gaussian landscape by introducing additional 
parameters and/or by scanning the "zero-mode wavefunction" 
of the Higgs boson. 
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are determined by integrating out the right-handed neu- 
trinos, giving the effective interactions 

{CX^I (0)) h t,hh, where = (A^ A^/ A„)^^. . (f f ) 

The field is a scalar that couples to the in the Majo- 
rana mass term because the Standard Model gauge sin- 
glets Ui are not necessarily gauge singlets in the higher- 
energy gauge group. In general there may be more than 
one scalar 0, and we take ip'^ to represent the effective 
wavefunction of the set of scalars that generate a mass 
term for the n^, after acquiring an effective vacuum ex- 
pectation value (0) . 

We assume the Majorana mass matrix A^^ of the right- 
handed neutrinos fii is generated randomly, via overlap 
integration in analogy to A" and the other Yukawa matri- 
ces. There is more than one possible origin of the right- 
handed neutrino mass interactions, and the overlap inte- 
gration for A^^ may not have the same phase functions as 
yu,d,e,n coming from the F™. Thus, we drop p(y) for A^^ 
in the numerical simulations of Section IlIII In fact, it is 
possible that neutrino masses are determined quite dif- 
ferently than this, however our assumptions capture the 
spirit of the Gaussian landscape, which replaces flavor 
symmetry charges with overlap integration of localized 
wavefunctions. 

In order to generate an ensemble of vacua and sta- 
tistical distributions of the flavor observables, we scan 
each of the randomly and independently over B. This 
ansatz borrows a hint from the fact that instanton center 
coordinates can be chosen freely in multi-instanton con- 
figurations, and fermion zero modes are localized around 
isolated instanton centers. It should be emphasized, how- 
ever, that this intuitive ansatz is far from justified from 
the top-down perspective outlined above — further un- 
derstanding is necessary to develop a rigorous connec- 
tion between microscopic theoretical formulations and 
this ansatz for generating statistics of flavor observables. 
Nevertheless, we consider it worthwhile to study whether 
such an ansatz generates phenomenologically acceptable 
distributions of flavor observables or not, before expend- 
ing too much effort on rigorous calculations. This sim- 
plified toy model of the string landscape is what we refer 
to as a Gaussian landscape. 

The Gaussian landscape above differs from that of 
Refs. [1, 0] in the inclusion of the kinetic overlap inte- 
grals, Eq. (O. These were ignored in the toy model of 
Refs. [1,13 because it was guessed that they would not 
significantly affect the essence of the results. We here 
find, however, that they have important implications for 
obtaining large lepton mixing in the Gaussian landscape. 

III. NUMERICAL SIMULATION OF THE 
GAUSSIAN LANDSCAPE 



for the iiT" and A° into predictions for the observable 
masses, mixing angles, and CP-violating phases. On the 
other hand, it is straightforward to perform such compu- 
tations numerically. Thus, we study the Gaussian land- 
scape using a numerical simulation, generating a large 
ensemble of sets of flavor parameters in order to repre- 
sent their landscape distributions."' As described in Sec- 
tion |TT1 we populate the landscape by taking the various 
peak positions y"' to scan randomly and independently 
over the internal manifold. Meanwhile, we take the other 
parameters — the microscopic gauge coupling g*, cutoff 
scale {M^,L)^, widths da/L, phases and phase func- 
tion p{y), the number of right-handed neutrinos Nn, and 
even the internal manifold B — as given. 

Some of the parameters that we take as given may in 
fact be uniquely determined by the microscopic theory, 
while others will scan in the full string landscape. Our 
strategy can be seen as focusing attention on a certain 
slice of the landscape. This creates a handle for under- 
standing the results, and meanwhile avoids the problem 
of guessing how these parameters are distributed among 
vacua in the landscape. Although the choice of values 
for these parameters is arbitrary, this will not be a prob- 
lem if the values we choose are not atypical in the full 
string landscape. For example, we choose du/L so as to 
match the observed level of hierarchy between up and 
top quark masses. It could be the microscopic theory 
uniquely determines this width, in which case our choice 
can ultimately be checked against this prediction, or that 
anthropic selection is important in determining the hier- 
archy, in which case our choice automatically includes 
this effect, or that du/L scans with a broad distribution 
that includes the value we choose, in which case the ob- 
served level of hierarchy is simply an accident. As is 
indicated below, for many changes in the values of these 
fixed parameters, the patterns of flavor remain intact. 

Let us now describe our choices for these parameters. 
First, we choose D = 2 extra dimensions, taking their ge- 
ometry to be that of the square torus T^, each dimension 
having length L. As is shown in some detail in Ref. 
the distributions of flavor observables in Gaussian land- 
scapes are broadly independent of the geometry of the 
manifold B, while the number of extra dimensions D 
primarily affects a single detail of these distributions, in- 
volving the distribution of the most massive particle in 
each generation. In the latter case, D ~ 1 stands out 
from the others, so we expect our choice of geometry to 
closely describe the results of any choice D > 2 on any 
smooth geometry. Because it is unclear how CP viola- 
tion should enter the Gaussian landscape — for example 
via the phase function p(y) in Eq. PU]) . via wavefunc- 
tion phases in Eq. ([7]) , or by some other mechanism — 
we set the issue aside for the moment and focus on the 



An analytic understanding of the Gaussian landscape 
is made difficult in part by the complicated set of op- 
erations necessary to convert the theoretical predictions 



We are grateful to Nathan Moore for helping to improve the 
efficiency of our numerical algorithm. 
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CP-conserving flavor observables, choosing p(y) = 1 and 
= for all a. We set g^s = g^/M^L = 0.5 to visu- 
ally match the observed b, and r Yukawa couplings; in- 
creasing or decreasing goff merely shifts the set of fermion 
mass distributions to larger or smaller values (note that 
QcS may determine the gauge couplings of the Standard 
Model; and therefore anthropic effects may be important 
in selecting its value). 

We choose da/L = d/L — 0.1 for all fields except the 
down-type quark singlets and the lepton doublets, for 
which dd/L = di/L = dg/L = 0.3. The choice d/L = 0.1 
is important for setting the observed level of mass hi- 
erarchy between the heaviest and lightest generations. 
The logarithmic spacing between typical masses in these 
generations is proportional to (L/d)^ Q; the level of hi- 
erarchy that we observe can be seen as an accident or 
as anthropically selected to generate appropriate atomic 
structure (see for example Rcfs. Meanwhile, in 

order to obtain large lepton mixing, the lepton doublet 
width dg must be larger. Hence we set de/L = 0.3. We 
choose dd/L = 0.3 simply because we consider it plausi- 
ble that these fields should have the same width, as they 
sit in the same SU(5) representation, 5 = {d,i}- As we 
shall explain, it would not hurt the fit to observation to 
give down-type quarks the same width as all of the other 
fields. 

Finally, we assume there are Nn = 20 heavy right- 
handed neutrinos. This number is chosen mainly for 
concreteness. As we shall describe, iV„ = 3 does not 
provide a terrible fit to existing observations for some 
choices of model parameters, including those described 
above. However, the fit to observation improves, and be- 
come more robust to changing the values of other model 
parameters, when Nn is increased. 

As we have already described, the peak positions y"'^ 
are taken to scan randomly and independently over the 
interval [0,L]. Although there are 22 flavor observables, 
for D > 2 the number of scanning parameters exceeds 
this, and the Gaussian landscape makes no "hard" pre- 
dictions. Instead, one obtains (correlated) probability 
distributions for the flavor observables. These are re- 
flected in the binned results of a numerical simulation, 
displayed in Fig. [1] For reference. Fig. [1] also gives ex- 
perimental values in brackets. These have been obtained 
by usin g re normalization group flow to run Particle Data 
Group mean values or constraints up to the Planck 
scale, where, for concreteness, we assume no particle con- 
tent beyond the Standard Model and ignore the effect of 
the right-handed neutrinos (alternatively, one could ex- 
tend the Gaussian landscape to cover low-energy super- 
symmetry) . We here emphasize that Fig. [1] displays the 
statistical distributions of flavor observables among vacua 
in the Gaussian landscape, and does not include cosmo- 
logical or anthropic selection effects, the latter of which 
we anticipate would modulate at least the distributions 

of XuM.e- 

The distributions in Fig. [T] are all quite broad. This is 
undesirable, from the point of view of experimental veri- 



fication, but may be the reality of fiavor in the landscape. 
Although the distributions are broad, the observed pat- 
terns among flavor parameters are present. For example, 
although the distributions of quark and charged lepton 
masses have considerable overlap within any given family, 
correlations among these masses are such that they are 
typically hierarchically separated. Furthermore, quark 
mixing angles are typically small, and Sfg^^ tends to be 
smaller than the others. It should be emphasized that 
the mass hierarchies and small mixing angles of the quark 
sector exhibit both the pairing structure and generation 
structure of the Standard Model, whereas neither of these 
is input directly into the Gaussian landscape. 

Note that the distribution of the d quark Yukawa cou- 
pling (and that of the electron) does not range over such 
small values as that of the u quark. This reflects a re- 
duced hierarchy, which is a consequence of the broader 
widths assigned to the down-type quarks (lepton dou- 
blets). If down- type quarks were assigned width d, their 
distributions would match those of up- type quarks. As 
they are, the distributions do not provide a bad fit to 
observation. However, if we consider an extension of the 
Gaussian landscape to include low-energy supersymme- 
try breaking, then tan/3 ^ 0(10) would seem to provide 
an improved statistical fit. Further investigation of this 
possibility is beyond the scope of this work. See also 
footnote m 

The mechanism that generates small mixing angles in 
the quark sector — which stems from the up-type and 
down-type Yukawa overlaps sharing the localized wave- 
functions of the quark doublets and Higgs boson — fails 
in the lepton sector, due to the broad widths of the lep- 
ton doublet wavefunctions. Ref. found that broad- 
width lepton doublets to be necessary but insufficient: 
without large CP-violating phases Va the relevant terms 
contributing to the PMNS matrix canceled. We find, 
however, that large CP-violating phases arc unnecessary 
if the kinetic terms K^j are determined by overlap inte- 
gration, as in Eq. The distributions of lepton mixing 
angles are broad, even with r = 0, so long as d^/L > 0.3. 

Ref. 3 found low-energy neutrinos to exhibit enor- 
mous mass hierarchy when Nn = 3. This hierarchy ulti- 
mately stems from rank reduction in the mass matrix C, 
which has the effect of "adding" the hierarchies among 
the eigenvalues of A" to those of A^-'^. Qualitatively, the 
low-energy neutrino mass hierarchy is related to the hier- 
archy among the three largest values of A^ . / Am; , where 
Xm and Am; denote respectively the eigenvalues of A" 
and A*^. When Nn = 3, the hierarchy among these 
terms is quite large. However, a consequence of choosing 
a large number of right-handed neutrinos, iV„ 3> 3, is 
their mass eigenvalues become densely packed. This is 
not hard to understand: the overall range of hierarchy 
is proportional to (L/c?„)^, and as more masses are in- 
serted in this range their typical separation decreases. In 
turn, the three largest values of A^ . / A Mi exhibit reduced 
hierarchy as 7V„ is increased. Indeed, we see in Fig. [1] 
that the hierarchy among low-energy neutrino masses is 
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FIG. 1: Distributions of flavor observables for the Gaussian landscape on T^. We set p{y) = 1, = for all a, g^s = 0.5, 
d/L = 0.1, d^/L = 0.3, and assume A^^ = 20 right-handed neutrinos. All logarithms are base ten, Xa denotes the Yukawa 
coupling determining the mass of particle a, and Ci is the ith eigenvalue of the low-energy neutrino mass matrix C"; other 
notations conform to those of the Particle Data Group. The scale (0) has been set so that the median of C2 is one. The 
numbers in brackets represent experimental values [ij], run up to the Planck scale. See the main text for details. 



not much larger than that among the quarks and charged 
leptons. 

Furthermore, it seems correlations tend to prevent the 
low-energy neutrino mass hierarchy from becoming too 
large in any given vacuum. As an example of this, in 
Fig. [T] we display the distribution of (the logarithm of) 
TO2/TO3, the ratio of the middle to the most massive low- 
energy neutrino. Wc sec that typically is not much 
more than 10-100 times larger than mj. Note that, in 
the limit of hierarchical low-energy neutrino masses, 

^Am|/AmL„ =^ m^/m^ (12) 

Hence in Fig. [1] we provide the observational constraint 
on the ratio of Am? of solar to atmospheric neutrino os- 
cillations. The Gaussian landscape provides a reasonable 
fit to this observation. 

Fig. [T] refers to a specific choice of right-handed neu- 
trino wavcfunction widths as well as that of the 
symmetry-breaking ficld(s) rf^. However, we have per- 



formed numerical simulations like that above but with 
dn/L = 0.3 or dcfi/L ~ 0.3, and have found that when 
D > 2 and Nn ^ 3 the distributions of fiavor observ- 
ables are qualitatively unchanged from those in Fig.[TJ It 
should be noted that choosing Nn = 3 increases the sen- 
sitivity of fiavor observables to the distribution of right- 
handed neutrino masses, so that mj/jrig becomes sensi- 
tive to dn- A larger low-energy neutrino mass hierarchy 
becomes more likely in this case, with the best fit to ob- 
servation coming from choosing narrow c?„ . 

Let us now discuss CP violation in the Gaussian land- 
scape. As has already been remarked, it is unclear how 
CP violation should be introduced into the Gaussian 
landscape in order to best reflect the phenomenology of 
string compactification. We have explored the two most 
obvious choices — using a F-matrix phase function p{y) 
or wavcfunction phases — and, remarkably, have found 
that both leave the distributions of CP-conserving ob- 
servables qualitatively intact. Furthermore, at least for 
the simple cases we have considered, both give similar 
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FIG. 2: Distributions of flavor observables for the Gaussian landscape on . Notation and parameters are as in Fig. [T] except 
r = 0.3. Distributions of flavor observables that are not displayed above are essentially unchanged from when r = 0. 



predictions for the distributions of CP-violating phases. 

For illustration, in Fig. [2] we display results of a numer- 
ical simulation exactly like that behind Fig. [U but using 

= r = 0.3 for all a (we study uniform Va = r merely for 
simplicity). Some CP-conserving flavor observables are 
not displayed in Fig. [2] because we could not detect any 
changes in their distributions relative to Fig. [TJ Those 
observables for which we could detect changes are dis- 
played in the top two rows of Fig. [51 It is seen that 
these changes are not very significant, considering that 
our universe provides only one data point to sample each 
distribution. The four CP-violating phases are displayed 
in the bottom row. Although the distribution of (JC'km 
is peaked at zero, the value we observe does not appear 
atypical. Both increasing and decreasing r from r = 0.3 
tends to sharpen the peak at 5'^^^ = 0, the effect being 
rather mild for increasing r. (Increasing r also tends to 
broaden the distributions of the quark mixing angles and 
of the mass eigenvalues Ad^.e;, while at the same time 
pushing their medians to smaller values.) 

We have also performed a numerical simulation using 
r = 0, p{y) = e^'^'^yi+^^^Z^ (for clarity we make explicit 
the two coordinates yi and y2 of the internal manifold 
T^). The resulting distributions are almost indistinguish- 
able from those obtained by choosing r = 0.3, p{y) = 1 
(Fig. [2]), except for the distributions of ai and a2, which 
are completely flat, and those of some CP-conserving ob- 
servables, which more closely resemble the distributions 
in Fig. m 



IV. THE AFS APPROXIMATION 

Although computations involving overlap integrations, 
diagonalizing mass matrices, etc., are most conveniently 
carried out numerically, it is possible to understand an- 
alytically the qualitative features of the distributions of 
various flavor observables in Gaussian landscapes. Such 
an understanding develops intuition for how these dis- 
tributions depend on the values of parameters that were 
taken as given in the numerical simulations of Sect ion [TTTl 
and thus allows us to anticipate the results of other 
choices of parameters without having to repeat the nu- 
merical simulation. Furthermore, as was seen in Ref. 0], 
such an analysis can reveal how the distributions of flavor 
observables are related to various features of the geome- 
try of compactification. Thus, we here extend the analy- 
sis of Ref. U by taking account of the kinetic mixing in 
the Kf; matrices. 

r I 

In Ref. [J], the present authors found that the dis- 
tributions of Yukawa couplings in Gaussian landscapes 
(using the basis of zero modes explored there) are almost 
the same as those predicted by models with approximate 
abelian flavor symmetries (AFS). In the AFS approach 
to flavor [12], one assumes each of the Standard-Model 
fermions has its own U(l) charge. If the ith anti-up quark 
has a U(l) charge ai > 0, the jth left-handed quark dou- 
blet has a charge bj > 0, and the Higgs boson is neutral 
under the U(l) symmetry, then their Yukawa coupling is 

Ar, =ff„e'^'e''^ (13) 
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where the small parameter e ^ 1 comes from a small 
breaking of the U(l) symmetry, and the gij are coeffi- 
cients of order unity, presumably determined by a higher- 
energy theory. Ref. Q showed that Yukawa couplings 
have the same structure in Gaussian landscapes, when 
the small suppression factors e"'' and e''^ are interpreted 
properly, even though there is no approximately pre- 
served U(l) flavor symmetry. This allows for a simple 
analytic understanding of quark and lepton masses and 
mixings. Below we show that a similar analytic under- 
standing of the predictions of Gaussian landscapes is pos- 
sible when including non-trivial kinetic mixing K°j and 
increasing the number of right-handed neutrinos. 

To proceed, we introduce a series of simplifying ap- 
proximations First, we ignore the fact that B is 
a compact manifold and instead treat it like Euclidian 
space when calculating overlap integrals. This ap- 
proximation becomes more accurate in the limit where 
wavefunctions are peaked nearer to each other — i.e. when 
the Yukawa couplings become large — but is surprisingly 
effective at describing the qualitative shapes of distribu- 
tions of Yukawa matrix elements, as was confirmed by a 
number of numerical simulations in Ref. Further- 
more, for simplicity we set p{y) — 1 and = (recall 
that the CP-conserving observables are approximately in- 
dependent of the Ta when Va ^ 0.3). 

With these approximations, the relevant overlap inte- 
grals have simple analytic expressions. For example, the 
up-type Yukawa matrix is given by 



exp 



3^2 



(14) 



where we have used translational invariance over the 
extra-dimensional space M.^ (this is true in as well) 
to set y'^ = 0. Note that for D > 1, the are vectors 
on the space of extra dimensions, and terms like y"'^y''^ 
correspond to scalar products of vectors. The last term 
in the exponent is "statistically neutral;" that is, it is 
linear in the random peak positions j/"* and , and 
therefore can be positive or negative, with a contribution 
that largely cancels in the statistical distribution. 
The up-type Yukawa matrix has the AFS structure, 

(15) 



where the small factors e" and correspond to the first 
two terms in the exponent of Xfj in Eq. p4)l . while the 
gij comes from the third term. Again, the gij are ran- 
dom numbers that may be larger or smaller than unity, 
but are statistically neutral. We can always relabel the 
generation index i of the three independent anti-up-type 
quarks Ui, so that e" < ^2 < ^3- Similarly, the gen- 
eration index for the three quark doublets can be rela- 
beled, so that el < e\ < £3- Note that in approximate 



It is interesting that curvature of B may introduce a statistical 
bias. We do not pursue this possibiUty here. 



AFS models, e is a fixed-value symmetry-breaking pa- 
rameter, and the U(l) charge assignments and bj are 
picked "by hand," whereas in Gaussian landscapes the 
"AFS factors" e^* < < and e\ < 4 < 4 follow 
statistical distributions. The quantity {L/dY sets the 
ordcr-of-magnitude of hierarchy in Gaussian landscapes, 
A '--^ e~(^/'^^ , playing a role analogous to the symmetry- 
breaking parameter e in approximate AFS models. 

For relevance to the results of Section IIIIl we as- 
sume the down-type quark wavefunctions have a rela- 
tively broad width dg ~ 0{L), while all other fields have 
the same (narrow) width d <^ L. The down-type Yukawa 
matrix then has the structure 



(16) 



The broad-width wavefunctions (pf are simply evaluated 
at the peak of the overlap, which occurs between the lo- 
calization of the Higgs boson y = y^ = and that of the 
quark-doublet wavefunctions y = y'^^ , and can be seen to 
introduce a set of random numbers of order the down- 
type wavefunction normalization. Thus, the down-type 
Yukawa matrix Af^ also follows the structure of approxi- 
mate AFS models. 

So far, our discussion has mirrored that of Refs. 0, 
However, that work ignored the effect of kinetic mix- 
ing. Applying the approximations described above, the 
kinetic mixing of quarks is given by 



if" 



exp 



4d2 



(17) 




1 











1 

















for a — u,d, q. Such a kinetic matrix can be put in canon- 
ical form by rotating quark flavors using a "triangular" 
matrix, for example the quark doublets are made canon- 
ical (with an appropriate normalization of ip"^) by 



(18) 



In the above, we have kept only leading order terms; note 
K'lj is of order (e^ e')^ for i j and thus the off-diagonal 
terms are small compared to unity. An analogous expres- 
sion applies to the up-type quarks. The down-type ki- 
netic matrices do not necessarily have small off-diagonal 
terms, yet they can be put in canonical form via a matrix 
of such "triangular" form. 

Note that these rotations do not affect the AFS struc- 
ture of the quark Yukawa couplings. The off-diagonal 
terms in the "triangular" rotation matrices for up-type 
quarks and quark doublets are negligible in this qualita- 
tive analysis, and so A" retains the structure of Eq. (|15p . 
Meanwhile, the "triangular" form of the down-type quark 
rotation essentially replaces the random terms (pfly^yij /2 
with linear combinations of such terms — this does not 
change the structure of A''. Thus, the quark sector of 
the Gaussian landscape retains its AFS structure. It is 
straightforward to see that the Yukawa eigenvalues are 
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roughly given by (i = 1, 2, 3), and (ej) 4 {j = 1, 2, 3), 
in the up-type and down-type sectors, respectively. 

As has been mentioned, the AFS factors e" and arc 
described by correlated statistical distributions reflect- 
ing the random scanning of Gaussian wavefunction peak 
positions. Their statistical distributions encapsulate the 
hierarchical suppression that must be inserted "by hand" 
using U(l) flavor charges in AFS models. Rcf. Q explains 
how the geometry of B affects the distributions of the 
AFS factors. That analysis carries over directly to the 
Gaussian landscapes considered here, so for brevity we 
do not repeat it. Instead, we simply note that distribu- 
tions of flavor observables are qualitatively independent 
the geometry of the internal manifold. 

Let us now turn to quark-sector electroweak mixing. 
Keeping in mind that kinetic mixing docs not change 
the structures of A" and A'' in Eqs. p5|) and (flBl) . the 
CKM mixing angles arc given by the difference between 
the diagonalization angles 9ij ^t/^'j < j) for tti6 up- 

typc Yukawa matrix and 9ij ^ (ff/Cj)^ < j) for the 
down-type Yukawa matrix. Qualitatively, the difference 
is dominated by the larger mixing angle, and hence 

^S™ - - A,./A,, ^ [\uj\u,)i ■ (19) 

The last two expressions relate the typical size of CKM 
mixing angles to typical ratios of quark mass eigenval- 
ues, due to all of these observables being qualitatively 
determined by the same AFS factors. 

The lepton sector can be analyzed similarly. Now the 
lepton doublet has broad width dg, while the other zero 
modes have narrow wavcf unctions. Thus, the charged- 
Icpton and neutrino Yukawa matrices have the structure 

- {4)-'^%=r^/2 (20) 

K, - ie7)iv%=y^./2. (21) 

These matrices have essentially the same structure as 
Xfj . Likewise, kinetic mixing in the lepton sector is either 
negligible (for electroweak singlets) or simply mixes the 
random factors cpj. Hence, the mass eigenvalues in the 
charged- lepton sector are qualitatively given by (ei)^, 
(e|)3, and (ef)*, where again we relabel generation in- 
dices so that el < e| < £3. Their distributions are the 
same as those of the down- type quarks. 

Because the right-handed neutrino mass matrix A^^ 
involves the overlap of two right-handed neutrinos, it is 
well approximated as a diagonal matrix, 

X^J ^ ef'S,, . (22) 

The diagonal elements ef^ may be small due to 
small right-handed neutrino overlap with the symmetry- 
breaking field (j). The seesaw mass matrix of left-handed 
neutrinos then takes the form 

a, - £ ^l\y=y-./2 ^%=y-./2 (^D V^f • (23) 
k=l 



The distribution of AFS factors - e-*?'"')'/^^' cov- 
ers a logarithmic range proportional to {L/dY . If, as we 
have assumed, the Majorana mass matrix A^ is deter- 
mined by overlap integration of right-handed neutrinos 
with a Gaussian wavefunction for (p, then the diagonal 
elements ef^ are also distributed over a logarithmic range 
proportional to (L/d)'^. The seesaw mass matrix Cij in- 
volves summation over all right-handed neutrinos, but its 
three eigenvalues are essentially determined by the con- 
tributions with the three largest values of Ck = (ffc) V^fe^- 
Although collectively the iV„ individual range over a 
logarithmic scale proportional to {L/d)'^, the three largest 
values strongly depend on Nn- As the number of right- 
handed neutrinos Nn is increased, typically the three 
largest values of Ck will become bunched near the largest 
possible value of Ck ■ In turn, the three eigenvalues of Cy 
experience reduced hierarchy as iV„ is increased. This 
is in good agreement with the measured ratio of Am^ 
coming from solar and atmospheric neutrino oscillations. 

Finally, we consider the distributions of the PMNS 
mixing angles. Because the different elements of a given 
row of A*^ and C are comparable to each other (due to 
broad- width lepton doublet wavefunctions), large-angle 
left-side rotations are necessary to diagonalize these ma- 
trices. Naively, these large-angle rotations would cor- 
respond to large PMNS mixing angles.^ However, the 
present authors concluded in Ref. (3] (where kinetic mix- 
ing was ignored) that broad-width lepton doublet wave- 
functions are insufficient to generate large PMNS mixing 
angles. The reason is, the broad-width wavefunctions 
tp^j [y) are not very different for different zero modes (dif- 
ferent j). As a result, the diagonalization angles for the 
left-handed charged leptons and those for the left-handed 
neutrinos are highly constrained, and an accidental can- 
cellation occurs in the PMNS matrix. Thus, it was con- 
cluded that some other ingredient (s) — such as large CP- 
violating phases — must be included in order for large lep- 
ton mixing to be typical in Gaussian landscapes. 

However, Ref. Q ignored kinetic mixing. As has been 
remarked, rotating leptons to a canonical basis of zero 
modes does not change the qualitative structure of the 
Yukawa matrices A*^ and A", nor the mass matrices A*^ 
and C . However, in such a basis the three independent 
zero modes of each family are orthogonal to one another, 
which is quite unlike the situation imagined in Ref. (3], 
where these zero modes were rather homogeneous and 
similar. The effect of this is to prevent the accidental 
cancellation that was observed in Ref. [1]. Thus, PMNS 
mixing angles obtain their "naive" values, comparable to 
ratios of different elements in a given row of A*^ or A", 



Note that this argument does not apply to the down-type quark 
mass matrix, with regard to the CKM mixing angles, because A'' 
has hierarehieal rows. The different elements of a given column 
of A"* are comparable, and this translates to large-angle right- 
side diagonalization rotations, but such rotations do not enter 
into the CKM matrix. 
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which arc typically order unity. Numerical simulations 
in Section [IIII confirm this idea. 

V. CONCLUSIONS 

The idea that characteristics of our universe are se- 
lected from a landscape of vacua has proved quite 
provocative. The major criticisms of this hypothesis can 
be classified into three groups: i) the landscape combined 
with eternal inflation implies a multiverse with diverging 
spacetime volume, expectations for observables depend 
on how this diverging volume is regulated, yet it is un- 
clear how this should be done; ii) expectations for ob- 
servables depend on anthropic selection effects, however 
it is unclear how to define our reference class, and at the 
same time the calculation of anthropic effects tends to 
be very complicated; and iii) our universe provides only 
one data point to compare to the predicted distribution 
of each scanning parameter, making landscape models 
seem difficult to falsify. 

On the other hand, the flavor physics of a landscape 
tends to avoid each of these criticisms. Although the cos- 
mology of eternal inflation and the spacetime measure of 
the multiverse (i) strongly influence the distributions of 
cosmological observables such as the observed primor- 



dial cosmic density contrast, insofar as the landscape is 
large enough that low-energy physics is decoupled from 
inflationary dynamics, flavor physics should be indepen- 
dent of these issues. Furthermore, although it is expected 
that anthropic selection (ii) plays an important role in 
determining some flavor parameters, such as the u, d, 
and e masses, most flavor parameters seem decoupled 
from conditions necessary to allow for the evolution of 
observers. Finally, although we have only one universe 
in which we can perform measurements (iii), there are 
almost two dozen (possibly correlated) flavor observables 
that should be described by the landscape model, allow- 
ing for a more discerning statistical analysis. Thus, we 
consider the issue of flavor to be ideal for applying the 
landscape hypothesis. 
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